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I. 


INTRODUCTION 


The corrugated horn has been established as an antenna with 
low side and back lobes, rotationally symmetric patterns (for square 
pyramidal and conical horn shapes), and broad-band performance[l -8] . 
These properties make this horn useful for many applications; including 
the one presently under study, i.e., as a radiometer antenna. A ray 
optics model of the dominant radiation mechanisms of conventional horn 
antennas is shown in Fig. la. In the corrugated horn (shown in profile 



CONVENTIONAL 
HORN GEOMETRY 

(a) 


CORRUGATED 
HORN GEOMETRY 

(b) 


Fig. 1. a) Conventional horn geometry b) Corrugated horn 
geometry 


in Fig. lb) the corrugated surface (with capacitive surface impedance in 
the E-plane of the horn) serves to reduce or eliminate the fields associ- 
ated with ray "a" of Fig. la. This in turn reduces or eliminates the usu- 
al high E-plane si delobes. The influence of the corrugation shape and 
density on the scattering from the onset of the corrugations and on the 
losses in the horn has been ignored in previous work. The parametric study 
reported here attempts to generate further understanding of the operating 



principles of the corrugated horn and to establish design criteria for 
the construction of practical corrugated horn geometries through a de- 
tailed analysis of the groundplane corrugated surface junction. The 
analysis includes studying the surface currents flowing on the corru- 
gations and the loss in the corrugations as well as the scattering from 
the onset of the corrugations. The usual corrugated horn requirements 
are for 8 or more corrugations per wavelength and a corrugation depth 
"d" between 0.25x and 0.5x.* In addition, it is usually specified that 
the tooth thickness ("t" in Fig. lb) should be much less than the corru- 
gation width ("W" in Fig. lb). In order to study the effects of varying 
these parameters, an integral equation solution was used to find the radi- 
ated fields and the surface currents associated with the groundplane cor- 
rugated surface models of Fig. 2. A brief description of the integral 
formulation is included in the next section and the results are presented 
in the following section. 

The groundplane-corrugated surface models snown in Fig. 2 were 
chosen for their similarity to one wall of a sectoral corrugated horn. 

In the corrugated horn, the illumination of the corruqations is by a cy- 
lindrical wave due to the diffraction from the horn-waveguide junction 
In the model used, the magnetic line source pair provides a similar cylin- 
drical wave illumination while at the same time allows placing a null in 
the far field in the direction of the adjacent edge of the groundplane. 

Then one need not match over the entire groundplane but only over the il- 
luminated part, thus saving the limited number of match points for the 
corrugated surface side of the model. The remaining match points are di- 
vided among the 20 corrugations. For all corrugation depths where cut-off 
operation is obtained (i.e., 0.25 i d/X < 0.5 for the square corrugation), 
the energy is forced away from the corrugated surface. Since the energy 
is forced off the corrugations, the corrugated surface matching points may 
be terminated without affecting the currents in the corrugations near the 
junction. Thus the surface model need not be closed. This conclusion was 
verified by finding the currents and scattered fields associated with the 
open model shown and a similar closed model and observing no significant 
differences. Therefore, this finite model is a good approximation to an in- 
finite groundplane corrugated surface junction as long as the corrugated 
surface is operated in the cut-off mode. 


*This depth would be in free space wavelengths for the parallel plate 
waveguide geometry of the corrugations. If the sides of the slot are 
terminated in a conducting plane, then these electrical lengths would 
correspond to those of the TE]g waveguide mode. Thus, some improve- 
ment in the operating bandwidth of the pyramidal corrugated horn could 
be achieved if the depth of the first few corrugations are appropriately 
tapered. 
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II. 


METHOD OF SOLUTION 


The surface currents on the corrugations and the scattering by the 
groundplane-corrugated surface junction of Fig. 2 were found using the H-field 
formulation for the transverse electric field case discussed b^ Harrington[9]. 
For this case, there is only a z component of magnetic field, H, and a tan- 
gential component of surface current, J(J = n x H, where n is the unit normal). 
At any point, the total magnetic field, H z , is the sum of the incident and 
scattered magnetic fields, H n z and H| respectively. The scattered field is re- 
lated to its source, the surface current, by 

(1) H| = U z • V x J JG d F 

where J = J dp = -[H z ] dp when H z is evaluated on C + (just outside the 

C _ 

contour C where the surface current, J, flows - the interior of C lies 
on the left side of <h') . 

and G is the two-dimensional Green's function GCp'.p' 1 ) = l/4jH^'k( Itr-p -1 1 ) . 

The resulting integral equation 

(2) J = -[H 7 + U 7 • v x f JG dP 1 ] 

22 J C + 

is solved for the surface current, J, by point matching using pulse basis 
functions. This integral equation reduces to the matrix equation 

(3) [ ] [fn] = [g m ] 

- - (2) 

where [e-mn] = j/4k AC n (n • R)H] (k|F m -Pn|) 

= matrix of coupling coefficients between the m^ and the n^ 1 
segments at p m and p n on the surface contour. 

R = — - — - = unit vector between n^* 1 and m^ points, 

I P"m - Fnl 

AC n = J n dr = current moment at the n t(l point, 

[f n ] = column vector of unknown surface current on n** 1 segment at 
the points (x n ,y n ) on the surface contour, 
anc * [9m] = column vector of incident field (= -H z ) at the points (X m ,Y m ) 
on the surface contour. 

This matrix equation has veen solved for the surfaces shown in Fig. 2 using 
a Crout[10] matrix inversion subroutine. Other more powerful inversion 
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subroutines which included pivoting and iterative improvement were used 
for comparison but provided no improvement beyond the 3 to 5 place accu- 
racy obtained from Crout's method. These other methods were not used 
regularly because they required nearly twice the storage since both the 
matrix and its inverse had to be stored. The subroutine listing for the 
Crout method is included in Appendix A along with the rest of the computer 
program. 




Fig. 2. Surface model used to find the scattered fields and surface 
currents associated with a groundpl ane-corrugated surface 
junction. 
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III. RESULTS FOR SQUARE CORRUGATIONS 

The parameters of the square corrugations of Fig. 2a which 
were varied in this study include the corrugation depth, the corruqation 
density and the corrugation shape. The notation used to relate these 
parameters to the .model and the range of values considered for these 
parameters are shown in Table 1. 

TABLE I 

PARAMETERS CONSIDERED IN STUDY OF SQUARE CORRUGATIONS 


PARAMETER 

NOTATION 
in Fig. 2a 

VALUES CONSIDERED 

Corrugation depth(in wave- 

d/A 

0.25, 0.3125, 0.375, 0.4375, 0.5 

lengths) 



Corrugation density (in 
corrugations per wavelength) 

NCOR = 1/TC 
or N/A 

4, 6, 8, 10, 12 

Corrugation shape! ratio of 



corrugation gap width to 

W/TC 

0.5, 0.6, 0.7, 0.75, 0.8, 0.9 

corruqation period) 

i 



Many combinations of these parameters were studied to ascertain the influence 
of each on properties such as the scattering from the groundplane-corrugated 
surface junction, the surface current flowing on the corrugations and also 
the power loss in the corrugations. 

Because of computer storage limitations, only 200 matching points 
could be used for each surface. These points were divided to allow for at 
least 8 matching points per wavelength along the surface. This matching 
point density is adequate if the fields in the corrugations are primarily 
those of the TEM mode. The verification of the dominant status of the TEM 
mode in the corrugations is contained in the higher order mode study in- 
cluded in Appendix B of this report. 

Of the 200 matching points available, ten were used on the ground- 
plane side of the model of Fig. 2a. The locations of these points are in- 
dicated in Fig. 3a. The physical optics current on an infinite groundplane 
is shown in amplitude and phase in parts b and c of Fig. 3. Also shown in 
Figs. 3b and 3c by the dashed lines are the amplitude and phase of the surface 
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current which exists at the midpoints of the surface segments of Fig. 3a. 
This demonstrates that the current model used in the calculations is an 
accurate model of the true surface current. 

The fields associated with the groundplane-corrugated surface 
junction are shown in Fig. 4 for corrugation depths from 0.25a to 0.4375a. 
The results shown are for a corrugation density of 8 corrugations per 
wavelength (i.e., corrugation period = TC = a/8) and for a corrugation 
shape ratio (ratio of corrugation gap width to corrugation period) of 
W/TC = 0.75. The patterns shown in Figs. 4a and 4b are the radiation 
zone total and scattered magnetic field patterns produced by the magnetic 
line sources of amplitude l/z 0 and the corrugations. The scattered field 
phase is shown in Fig. 4c. Notice that the direction of maximum radi- 
ation intensity varies from about 65° above the corrugated surface for 
the 0.25a depth to about 35° for the 0.4375a depth. All of the other 
corrugation densities and corrugation profiles showed essentially the 
same patterns at the same corrugation depths and are therefore not in- 
cluded. The smoothly varying phase pattern indicates that the phase 
center of the scattered field is nearly at the origin of the coordinate 
system (the junction between the groundplane and the corrugations). 

This range of angles for this scattered field maximum explains 
the peculiar behavior of the small corrugated horn discussed previouslyDU. 
This horn had a pattern which showed a slight frequency dependence. How- 
ever, as the end of this horn was cut-off, the pattern of the remaining 
corrugated horn became quite frequency dependent. The geometry is shown 
in Fig. lb. For the original corrugated horn, the scattered field maxi- 
mum illuminated the opposite wall of the horn and weakly illuminated 
the aperture edge (the aperture edge was at about 43°). However, after 
the modification, the edge of the horn occurred at an angle of 50° and thus 
was more strongly illuminated. The fields diffracted by this edge 
phased destructively and constructively as a function of frequency with the 
desired radiated field and thus caused the frequency dependent patterns. 

The magnetic field intensity associated with the rays diffracted by the 
corrugated surface junction may be expressed in the form 


n i o"J’ kr 

(5) H° = H 1 2-y- F(e) 

where H 1 is the magnetic field incident on the junction, 
r is the range as measured from the junction, 
e is the angle shown in Fig. 2, and 
F(e) is the complex pattern factor shown in Figs. 4b and 4c. 

It is noted that the phase is a slowly varying function of 0 and thus pre- 
sents no difficulty for the application of GTD. This information is now 
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sufficient to compute the diffracted fields associated with ray "b" of 
Fig. lb. The far field of the ray in the direction of the horn edge is 
readily obtained from the preceding equation. Next the surface discon- 
tinuity is replaced by a line current whose fields have that magnitude. 

This line current then illuminates the 90° corner at the end of the horn 
wall and the diffracted fields are now computed using the techniques of 
the Geometrical Theory of Diffraction. If the fields in the deep shadow 
are to be computed, then the diffraction from the thick wall is computed 
by representing it as a pair of 90°wedges[l 1 ]. While these computations 
have not been carried to completion at this time, they represent a straight- 
forward implementation of accepted GTD practices and will yield accurate 
results. The major goal at this time is to establish the condition for 
which this edge diffraction is negligible and this would occur when the 
angle to the opposite edge is less than 20°-25°. 

As mentioned earlier, another property of the corrugated surface 
which is of interest is the rate of decay of the surface current flowing 
on the corrugation walls. The decay in the amplitude of the surface cur- 
rent is due to the energy being forced away from the corrugations ana not 
caused by power loss in the conductor. The loss will be diccussed later. 
Figure 5 shows the normalized surface current versus corrugation number 
for a surface with 8 corrugations per wavelength and a profile ratio of 
W/TC = 0.75 at several corrugation depths. The surface current plotted is 
the current which exists at the bottom of each of the 20 corrugations (point 
"B" in the insert) normalized with respect to the surface current which 
exists at the same X-coordinate on an infinite groundplane with the same 
sources acting (point "A" in the insert). This normalization removes the 
range dependence of the incident cylindrical fields and also selects the 
maximum current which exists on the surface of the corrugations. 

Since to a good approximation only a TEM mode exists in the cor- 
rugation, the currents on the teeth at any point are given by 

(6) |J y | = J s cose (d + y) 

where y assumes negative values and 3 = 2u/a 0 is the propagation constant 
of a TEM wave in the parallel plate waveguide region. This has also been 
established by all of the computations that have been made. 

Notice that the most rapid decay is obtained for the 0.25 a case and 
as expected no decay occurs for the 0.5a case. The variations of the current 
near the end of the surface are caused by fields reflected from the termi- 
nation of the structure and may be ignored. This property (the decay of the 
surface current) has been examined for other corrugation densities and pro- 
files. Figures 6 and 7 show the surface currents (normalized as above) exist- 
ing on 0.25x and 0.375a deep corrugations for various corrugation densities. 
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Radiation pattern associated with a groundplane-corrugated surface 
junction illuminated by a cylindrical wave: a) total fields, 
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Decay of surface current on a corrugated surface due to the 
energy being forced away from the corrugations. 
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Decay of surface currents on 0.25 a deep corrugations as a function 
of corrugation density. 




Decay of surface currents on 0.375A deep corrugations as a function 
of corrugation density. 
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Decay of surface currents on 0.25X deep corrugations as a function 
of corrugation shape. 



Decay of surface currents on 0.375 a deep corrugations as a function of 
corrugation shape. 
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Power loss in a corrugated surface for various corrugation depths. 
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Power loss in 0.25x deep corrugations from the origin to "X" (normalized 
with respect to the loss in a finite segment of an infinite groundplane 
from the origin to "x") for various corrugation densities. 
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Power loss in o.375\ deep corrugations from the origin to "X" (normalized 
with respect to the loss in a finite segment of an infinte groundplane 
from the origin to "x") for various corrugation densities. 




These results are plotted as a function of the X-coordinate measured in "cm" 
from the onset of the corrugations (the X and 2x points are also shown on 
this scale, A = 8 cm). Corrugation densities ranged from 4 to 12 corrugations 
per wavelength and the results for each corruqation depth showed essentially 
no dependence on the corrugation density. The dependence of the surface cor- 
rent on the corrugation profile was also investigated. The corrugation pro- 
file was specified by the ratio of the corrugation qap width (W) to the cor- 
rugation period (TC) and the corrugation density was 8 corrugations per wave- 
length. The results are shown in Figs. 8 and 9 for 0.25 a and 0.375 A deep 
corrugations, respectively. The range of W/TC ratios considered was from 0.5 
(thick metal vane between corrugations) to 0.9 (very thin metal vane). 'Beyond 
approximately 1/2 a from the onset of the corrugations ( i . e . , x > 4 cm for the 
A = 8 cm case shown), the decay of the surface current per wavelength is near- 
ly constant (■»- 6 dB/A at d/A = 0.25 and ■». 4 dB/A at d/A = 0.375). Thus in a 
practical situation, one could use very thin vanes near the onset of the cor- 
rugations and thicker vanes (which are easier to construct) further from the 
onset of the corrugations . These curves then establish the length of corru- 
gated horn required to reduce the edge diffracted fields. Since the current 
density at the top of the teeth is almost zero (see Appendix B) only a few 
teeth would be required here. For the 3/8A depth, the currents at the top of 
the tooth would be 3 dB below those at the bottom. Thus the illumination of 
the distant end of the horn would be -15 dB for x T 16 cm = 2A for W/TC > 0.7. 
Thus the fields of the diffracted ray ("a" of Fig. 1) that phase constructively 
and destructively with the geometrical optics fields are reduced to a negligible 
value with a surface that is two wavelenqths in extent. 

Another property of the corrugated surface which was investigated is 
the power loss in the corrugations (both the total loss and also the location 
of regions where the maximum loss occurs). In an effort to make the results 
more general, the corrugated surface loss is normalized with respect to the 
loss in an equal size segment of an infinite qroundplane of the same material. 
This eliminates the need to assume a surface resistivity for the material. The 
results shown in Fig. 10 show the normalized loss in dB versus corrugation num- 
ber "N" for various corrugation depths. The loss in the corruqated surface is 
the sum of the loss in each corrugation from the onset of the corruqations fat 
the origin of the coordinate system) up to and including the loss in the W™ cor- 
rugation. The loss in the groundplane (used for normal i zati on) is the loss in 
the region from the origin up to the X-coordinate of the end of the N^ 1 corru- 
gation. Notice that most of the loss occurs in the first few corruqations and 
also that it appears that all of the curves (except the 0.5a deep case) would 
cross the 0 dB line (equal loss in corrugations and groundplane) after a reason- 
able nurrber of corrugations (the 0.3125a and 0.375a deep cases cross in fewer 
than 20 corrugations). The power loss is also dependent upon the corrugation 
density and shape. Figures 11 and 12 show the power loss in 0.25a and 0.375a 
deep corrugations for various corrugation densities. This loss is the loss in 
the corrugations from the origin to "X" normalized with respect to the loss in 
the same length segment of an infinite groundplane. Notice that the lower den- 
sities are less lossy than the high density surfaces. This result is reasonable 
in view of the fact that the surface current decay per wavelength is nearly inde- 
pendent of the corrugation density (c.f., Figs. 6 and 7). Thus the higher corru- 
gation density with its greater surface area should be more lossy. The influence 
of the corrugation shape on the power loss may be seen in Figure 13. 
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Shown here is the loss in 20 corrugations (from x = 0 to x = 20 cm) nor- 
malized with respect to a 20 cm segment of an infinite groundplane ver- 
sus the corrugation depth (D/A) for various W/TC ratios. The W/TC = 0.5 
case corresponds to a thick metal vane between corrugations while the 
W/TC = 0.9 case corresponds to a very thin vane. Notice that the thin 
vanes have lower loss than the thick vanes and that there exists a range 
of depths over which the loss is minimized. 


PLANAR SURFACE 
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Fig. 13. Power loss in 20 corrugations versus corrugation 
depth for various corrugation profiles. 
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The general conclusion which might be drawn at this point is that 
corrugation density and tooth shape can be selected to optimize the 
performance of a corrugated horn by a) choosing as few corrugations as 
four/wavelength, b) by making the vanes as thin as is practical. However, 
there is one final parameter that must be considered. The impedance of 
the horn would be related to the magnitude of the wave reflected back 
toward the source. The computer programs developed here are not well 
designed to evaluate this reflection coefficient. However, they can be 
used to obtain an indication of this parameter. The current on the 
groundplane region of Fig. 2a can be expressed in the following form: 

(7) J = J s + J 1 + J 2 

where 7 is the current density of the source when placed on an 
_ s infinite conducting plane, 

J-] is the current density associated with the wave reflected 
_ by the onset of the corrugations, and 

Jp is the current density associated with the waves reflected 
c by the end of the groundplane. 

Since the two line sources are phased to give zero far field in the 
e = it direction, and since the diffraction coefficient associated with 
the back scattered field at this edge is small, J2 may be neglected. 

When the observation position is sufficiently removed from the edge, 
the current density 3] will decay as 1//F since it is caused by the 
radiation of a line source. Thus the equation for the surface current 
has the form: 

(8) J = I s + 2n x H q //F 

where r is the distance on the groundplane from the onset of corrugations. 
Since J 5 is known exactly, J] can be found with reasonable accuracy. Its 
representation as 2H 0 is, of course, only approximate. A plot of |H 0 | 
versus the corrugation depth is shown in Fig. 14 for the W/TC = 0.75, 

8 corrugations per wavelength case. As previously discussed, the a/ 4 
depth has a large reflected wave associated with it. The reflection 
from the junction also depends on the tooth shape and density. Figures 
15 and 16 show |H 0 | versus corrugation density and profile ratio re- 
spectively. Notice that at both x/4 and 3x/8 depths, the higher cor- 
rugation density and thick teeth (W/TC =0.5) have lower reflections 
associated with them. The W/TC = 0.9 points in Fig. 16 are shown with 
a dashed line as some difficulties were noticed for cases where matching 
segments were very close together. This is probably the case for the 
W/TC = 0.9 data. 

A general observation is that the impedance as a function of cor- 
rugation density would require as many corrugations as possible in order 
that the VSWR be maintained low over the largest possible bandwidth. 

Now we arrive at the optimum design of the corrugated surface in 
terms of expense, and operation. There should be from 8 to 12 corrugations 
over the first wavelength. These teeth should be as thin as is possible for 
low loss. Then the tooth density should be tapered to 4 per wavelength and 
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the tooth width should be increased to W/TC = 0.5. In addition, the first 
few teeth should be plated to further reduce losses. Also, for low loss, 
the waveguide and throat region of the horn should be plated regardless of 
whether the horn is corrugated or not. 

However, if loss is not an important consideration, then the reflec- 
tion coefficient and consequently the horn VSWR should be minimized by in- 
creasing the corrugation density to at least 8 corrugations per wavelength. 
The reflection coefficient can also be reduced by initiating the corrugation 
at a point further removed from the horn waveguide junction. However, one 
must be certain to initiate the corrugations before a second mode could be 
excited in the horn[4]. 



CORRUGATION DEPTH (d) (WAVELENGTHS) 

Fig. 14. | H 0 1 versus corrugation depth. 
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IV. 


RESULTS FOR VEE CORRUGATIONS 


The study of the vee corrugations shown in Fig. 2b was undertaken 
because there are situations where the square shape of the previous section 
is impractical or unsuitable. The thin vane of the square corrugation is 
impractical at frequencies much above 12 GHz because of the difficulty en- 
countered in machining or extruding the corrugations. The difficulty should 
be eliminated with the V-shape corrugation. Also, the square shape is not 
as suitable as the V-shape for an application where the surface needs to 
fold or collapse in order to take up less space. Using an asymmetric V- 
shape, an unfurlable corrugated horn antenna which would fold like a camera 
bellows is practical. The corrugated surfaces considered in this section 
have symmetric V-shape corrugations defined by the corrugation depth, d/A, 
and the corrugation density, NCOR (in corrugations per wavelength ) . The 
range of parameters presented for the V-shape is not as complete as for the 
square shape because of some problems encountered in calculating the surface 
currents for the deep corrugations and for high corrugation densities (NCOR 
8). These cases will not be reported here. It is believed that for these 
cases, the coupling coefficient between adjacent segments at the top or bot- 
tom of each corrugation is erroneous because of the very small distance be- 
tween the matching points. A sufficient number of cases are available to in- 
dicate that the V-shape does operate as a cut-off corrugated surface for cor- 
rugation depths from 0.3125 a to about 0.625A. Cut-off operation of a corru- 
gated surface requires a capacitive surface impedance. That the surface im- 
pedance of the vee corrugation is capacitive over most of this ranqe of depths 
has been verified using Harrington's results for the fields in a wedge shape 
waveguide with radial propagation[12] . The impedance looking into the wedge 
waveguide is: 

(9) Z = jZ 0 J 1 (kp)/J 0 (kp) 

where p is the radial distance from the bottom of the corrugation. For the 
deep, narrow corrugations being considered, this yields an operating range of 
approximately 0.382 .1 d/A £ 0.605. This operating band is narrower than the 
computer calculations indicated but then the wedge radial waveguide is only 
an approximate model for the V-shape corrugations. The results of the com- 
puter calculations are discussed next. 

Figure 17 shows: a) the groundplane matching points, b) the sur- 

face current magnitude, and c) the surface current phase in the source 
region. Figure 17b shows the physical optics current (on an infinite ground- 
plane) and the approximate current representation for the computer selected 
match points. Both the amlitude and phase are accurately modeled for this 
point selection. 

The properties of the vee corrugations which were studied include the 
radiation patterns of the surface model, the surface currents in the corru- 
gations, and the power loss in the corrugations. Figure 18 shows the radi- 
ation pattern (total magnetic field) of the surface model shown in the insert 
and in Fig. 2b for several vee corrugation depths. As with the square cor- 
rugations, the radiation pattern is nearly independent of the corrugation 
density (i.e., the pattern at a given depth is nearly identical for the 4, 6, 
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Fig. 17. Surface currents near the line sources used to model 

the groundp lane -corrugated surface junction a) matching 
points on groundplane, b) surface current magnitude, and 
c) surface current phase. 
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or 8 corrugations per wavelength cases) and the direction of maximum in- 
tensity decreases with increasing corrugation depth. 




Fig. 18. Radiation pattern of a groundplane-vee corrugation 
junction illuminated by a cylindrical wave. 


Other similarities between the V-shape and the square corrugations 
are found in the surface current decay and power loss. Figures 19, 20, and 
21 show the surface current (in dB) which exists near the bottom of the vee 
corrugations normalized with respect to the surface current which exists at 
the same "X" coordinate on an infinite groundplane with the same sources act- 
ing. Figure 19 shows this current versus the distance (in wavelengths and in 
cm) from the onset of the corrugations for two corrugation depths near the 
low end of the operating band. At the lowest frequency in the cut-off band 
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onset of the corrugations. 
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Relative surface current on 0.375 a deep vee corrugations versus the distance 
from the onset of the corrugations for several corrugation densities. 











Relative surface current on 0.5 deep vee corrugations versus the distance fror 
the onset of the corrugations for several corrugation densities. 
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Relative power loss in vee corrugations from origin to "X" (normalized with respect 
to the loss in 0-0( segment of infinite groundplane) for two corrugation depths. 



Relative power loss in 0.375 deep vee corrugations from origin to "X" (normal- 
ized with respect to the loss in 0+X segment of infinite groundplane) for three 
densi ties . 



(d/A = 0.3125), the surface current decay along the vee corruqated surface 
is nearly identical to that of the square corruqations at the lowest frequency 
in its band'(c.f., d/A = 0.25 in Fig. 5). As was the case for the square cor- 
rugations, the surface current decay is most rapid at the low end of the fre- 
quency band. Figures 20 and 21 show similar current curves for various cor- 
rugation densities at depths of 0.375a and 0.5A respectively. Notice that near 
the lower end of the operating band the surface current decay is independent of 
the number of corrugations per wavelength while near the upper end of the band 
there is some dependence on the corrugation density. 

Figures 22 and 23 show the relative power loss in the vee corrugations. 

The loss which is plotted is the loss in the corrugations from the origin to 
the "X" coordinate normalized with respect to the loss in the same length seg- 
ment of an infinite groundplane. Figure 22 shows the loss at the two lowest 
frequencies for the 8 corrugations per wavelength case. The loss here is higher 
than the square corrugations because of the higher surface current near the point 
of the vee and also because of the greater arc length in the V-shape surface. 

Figure 23 shows the loss at d/A = 0.375 for 4, 6, and 8 vee corrugations per wave- 
length. As with the square corrugations, the lower corruqations density has the 
lower loss. Because of the higher current (and consequently higher loss), if 
one is to use the V-shape corrugations in a situation where low loss is important 
something such as polishing or silver plating the first few corrugations should 
be considered. For an unfurlable horn, the first part of a horn might be rigid 
and use square corrugations to obtain lower loss while the remainder of the horn 
could be an unfurlable geometry using triangular teeth. Also, because the cur- 
rent in the corrugations become very small after 4 or 5 corrugations the loss in 
the corruqations will eventually become smaller than the loss in the groundplane. 

V. CONCLUSIONS 

A parametric study of square and vee corruqations has shown that both 
shapes operate as cut-off corrugated surfaces over subtantial bandwidths with 
relatively low loss. The study has also shown that the rate of decay of sur- 
face current is nearly independent of the corrugation density (the number of 
corrugations per wavelength and that the lower corrugation densities (4 or 6 
corrugations per wavelength) are desirable because of the lower loss. A sig- 
nificant part of the loss is confined to the first few corrugations and should 
allow one to treat this region with special care such as silver plating or special 
polishing of the surface when using corrugated surfaces in situations where low 
loss is important. In any event, the loss in the walls for most practical cor- 
rugated horns would be as low as (if not lower than) the loss in the same shape 
and size conducting wall of the same material with the same surface finish. The 
results also indicate that the VSWR of the horn will depend on the corrugation 
shape. For the groundplane-corrugated surface junction considered, the VSWR de- 
creased with increasing corrugation density and increased with thinner teeth be- 
tween the corrugations. Since the loss increases with increased corrugation den- 
sity, some compromise is required. Six to eight corruqations per wavelength near 
the onset of the corrugations and then decreasing the density to two to four per 
wavelength should be adequate. Again, since the loss was lower and the VSWR high- 
er for the very thin teeth,; some compromise is also indicated A corrugation width 
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to period ratio of W/TC = 0.75 should be a good compromise value if one v/ants 
to avoid the complexity of changing the corrugation shape along the corrugated 
wal 1 . 


While this work has yielded considerable insight about corrugated sur- 
face design, there remain several sets of important design data that can be 
obtained with a similar approach. For example, the sources on the qroundplane 
should be moved off the groundplane and the corrugated surface performance an- 
alyzed for various source locations. This would simulate the fields incident 
on the corrugations due to the diffraction from the horn-wavegui de junction and 
the onset of the corrugations on the opposite wall of the horn. Also, the ground- 
plane should include a corner similar to the horn-waveguide junction. Ideally, 
the entire horn should be considered. This has been programmed usinq the inte- 
gral equation approach but is of limited value because of the extremely small 
horn which can be handled due to computer stnrage limitations. Using image tech- 
niques and other procedures currently being developed at the ElectroScience Lab- 
oratory, it should be possible to consider 10 or more corruagtions in this model 
for VSWR computations . This number of corrugations should provide adequate re- 
duction of the surface current at the top of the last corrugation to obtain valid 
calculations for realistic horn throat geometries. For the corrugated horn pat- 
tern, a diffraction coefficient could be associated with the scattering from the 
onset of the corrugations using the integral equation approach. Then, using the 
techniques of the Geometrical Theory of Diffraction, one could proceed to obtain 
full horn patterns including the back lobe radiation. 
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APPENDIX A 


COMPUTER PROGRAM USED FOR CORRUGATED SURFACE ANALYSIS 

The main computer program and the subroutines required to solve the 
set of equations obtained from the integral equation of section II by point 
matching will now be briefly discussed. The program and subroutines are 
written in the Fortran IV language. The main program listing is at the end 
of this section. The input parameters to the main program are contained on 
two cards: the first specifies the corrugation shape, and the second the 

groundplane and corrugated surface dimensions. The data on the first card 
is grouped into five groups of ( F8 .3,12) fields. Each group locates a cor- 
ner (in F8.3 field) and specifies the number of segments (in 12 field) to be 
placed on the line from the previous point to the current point. Figure A1 
shows a typical period or cell from a corrugated surface of period TC = 1/ 
NCOR where NCOR is the corrugation density (in corrugations per wavelength ) . 


H - 0 


N5 SEGMENTS 

X = 0 XIP X4P I X5P - TC 


I I 111 


N I SEGMENTS 
N2 SEGMENTS 

H = HC 



N4 

SEGMENTS 


Lj 


N 3 SEGMENTS 


X = 2 TC 

I 


Fig. A1 


The XIP to X5P are expressed as ratios normalized to the corrugation period 
TC with X5P = TC. The heights of the comers are assumed to be either H = 0 
or H = -HC as shown. Additional information on these parameters is contain- 
ed in the discussion of the VSURDV subroutine. The remaining input parameters 
are shown in Fig. A2 and listed below: 

EPL locates the end point on the l_eft or groundplane side in cm., 

NGP specifies the number of matching segments to be placed on the 
groundplane side of the origin. 
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EPR locates the end point on the riqht or corrugated surface side 
of the origin (in cm.), 

WE is the electrical wavelength (in cm.), 

HCWE is the height of the corrugations in wavelengths, 

NCOR is the number of corrugations per wavelength, and 
KSET is an identifier for the set of data and is not used in anv 
cal cul ations . 


EPL (cm) x=0 WE EPR (cm) 



Fi g. A2 ' 


Also shown in Fig. A2 as a broken line is the path used for closing the body 
to check the validity of the open surface model. The body is closed by set- 
ting the last calling parameter in the CALL VSURDV (line 36 of the main pro- 
gram listing) to 1 rather than 0 (open body). The match points on the ground- 
plane and corrugated surface are determined by calling the GDPLDV and VSURDV 
subroutines in this order (lines 35 and 36 in the listing). These subrou- 
tines are discussed later. The output includes a list of the "X" and "H" co- 
ordinates of the match points, the direction of the surface normals, and the 
length of each segment as determined by these subroutines (line 39). Also 
listed as output is the incident field on each segment (line 74) and after 
the matrix inversion, the surface current on each segment (line 82). Also 
plotted are the magnitude (line 88) and phase (line 93) of the surface current 
versus segment number and the radiation zone magnetic field pattern versus the 
observation angle (line 122). In addition to these plots the magnitude and 
phase of the scattered field and the total field are listed versus e(line 12 ). 
Finally, the relative power loss is summed and listed versus the segment number 
(line 136). The Fortran listing follows. 
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1 r T* CrtSf GAUSSIAN I NTE6R A 1 T OH USFO TO Fill IN MATRIX 

2 C T | CAGf I.AlisMAfl iNTfGRAnON US| 0 TO Fill IM MATRIX 

•»> fOHI’U K lin<( .mr ,M(!Ol<Nl 

4 f ClMI'l f v rA/. f ro^MN,s^T 

f, CUMPU > T ,( H SI .n,M*HAA/P»nJC*MS 

6 cu^mon /po&/ PJ r 

7 DIMENSION Aprs ( 36 0 ) , Y ( 1 U ) ,HTOT ( 360 ) 

o rO M MON /MAW J GO/ nZl,UZ2 

g rO'iMoN/MA JM>V/ X(2^0),H(?30)»TN(?3Q)*SEGL<230> *OEL*Hc*G*LPL*EPR*WE 

in complex F INC <2^0 > «C (230,230) ,r (230 ) 

11 NDTM=230 

IP C V'E IS THE ELECTRICAL, WAVELENGTH 

13 C SOHRrt. LOCATIONS 
1M PZl=-4.0 

15 n/2--6.0 

16 RE AD <5,1 ) XlP,Ml,X2P«N2,X3P,M3,X4P*N4*X5PtN5 

17 1 PORN AT (5(F8.3,T2) ) 

16 RLAD< tin?) EPL OIGP. f_PR,WF,HCWF« NCOR*KSET 

19 2 FORMAT (FlO. 4, I*«3P10.4t 15*25/ , 15) 

20 WRITE (6,9997 ) EPL • EPR , wE t HC wF t NCOR , NGP 

21 9997 rORMAT ( * GROUND PL A NE FROM *,c 1 u# 4, * CM TO ORIGIN, VEE CORRUGATION 

2? PR FROM ORIGIN TO *,Flo,4,* C M, /» *jAVELENGTH = • , FlO . 4 , • CM*/* CORRUG 

23 3ATI0M HFI(ihT=» ,Fin.4, *WAVLLf NGTHS* /IS, ’CORRUGATIONS 

24 4PlR WAVELf M»TH f /I5* ♦ POINTS ON GROUNDPL ANE * ) 

25 PJ=3, 14159 

26 Pl2=PI/2. 

27 R0FG=160,/PI 

28 TC=WF/NCOR 

29 y 1 -X 1P*TC 

30 X2=X?P#TC 

31 y3=X3P*TC 

32 X4 -X4P* TC 

33 y5=X5P*TC 

34 MC=HCWE*WE 

35 CALL GPPLOV(NGP*MP ) 

36 CALL VSURQV(Xl,X2,X3,X4,TC*NltN2,N3,N4,N5,MP,0 ) 

37 MsMP -1 

38 U)RT TF ( 6 , 20 ) 

39 wRlTf (6,*.1)(M,Y(M),H(M),TN(M) ,SEGL( M) *M = f *N) 

4 0 20 FORMAT (T4, *M* ,T15 , • X* , T30* *H f * T45* • TN f ,T60* 'SEGL* ) 

41 21 MlRMftT (I5»4fl5 # 4) 

4? TF(N.LF.NOIW) GO TO 19 

43 WK I TF ( 6 , 1 6 ) N 

44 18 PORMATC N= * , I 5 , * EXCEEDS DIMENSION OF C») 

45 CALL EXIT 

46 19 CONTINUE 

47 G ^6. 2831853 /WE 

4ft STS= -CMPLX < 0,707107,0,707107 >/< 2.0 *SQRT( We > > 

49 njC=CMPLX ( 0,0,1 ,0 ) *G/4.0 

50 NMO=N-l 

51 Nr,3rN-3 

52 PPIF=0, 7853982 

53 C MATRIX FILL IN 

54 00 3661 IR=1,N 

55 no 3661 IC=1*N 

56 3661 C(IR,IC)=CO< IR,IC) 

57 r this completes t h e fillin of the matrix 

58 WRITE (6,1 22? ) N*WE 

59 1222 FORMAT ( 3H N=,I3*4H WE=,E15.8) 

60 c THIS FINDS THf XNclOENT F I ELp IOM THE N JTH SfGMENT 

61 00 455 Nj=ri , N 

62 xG = 0.5O*(x(NJ)+x( NJ + 1 ) ) 

63 r THE SIGN Ofi IMF INCIDENT FIELD HaS BEEN aOJUSTEO TO aGREE WITH 

64 C THE INTEGRAL EQUATION 

65 M6 = 0,5tl0» (H(NJ) ♦H(NJ*1) ) 


35 



66 


RHn = SQRT ( XG*XG+N6*H6> 


67 


PHI = A1AN2CH(m XG)*RDEG 


66 


rALL W1LLUM(RH0*PHI ,HINC,WE) 


69 


FINC(NJ)=-HIMC 


70 


F<NJ)=FINC(NJ> 


71 

455 

CONTINUE 


72 


WHITE (6,2947) 


73 

2947 

F0RMAT(/5H m * 1 0 X * 1 6H INCIDENT FIELDS ) 


74 


WRITE (6,2948) ( N J « F INC ( N J ) «NJ=lf N) 


75 

2946 

FORMAT (I5,E15.8«4H + J , E 1 5 , 6 ) 


76 


CALL CROUTN ( C t F « N « NDIM ) 


77 


WRlTr(7t9300) (u'SET,M,X(M) ,H(M) ,TN(M) *SEGL<M) ,F(N) fM=l»N) 

76 

9300 

FORMAT (2I5«4Ein,4,2El5«8) 


79 


nO 554 IKtlR = l,M 


60 


A AFrC ARS ( F (IKIJR) ) 


61 


ANF = 57 #296* A TAM? ( A JMAG ( F ( I KUP ) ) « REAL ( F ( I KUR ) ) ) 


6P 

554 

WKTTE( 6,5 53) I KUR i A AF » ANF 


63 

553 

FORMAT (• • , «F ( • , I4 t f ) = ’ • E15.0 t * AT ANGLES • » E 1 5 . 8 ) 

64 


HO 9553 IRR0=1 • N 


85 


INn=IRRO-l 


86 


Y ( 1 ) =C ABS ( F ( IRRO) ) 


87 


XRRO=FLOAT ( IRRO) 


66 

9553 

CALL PLOT(XPRO« Y*1 t JNDt5, 00,0.0) 


89 


00 9554 IRB0=1.N 


90 


IN0=IRR0-1 


91 


Y( 1)=57,255C*ATAN2( AIMAG(F ( IRRO) I .REAL (F ( IRRO) ) ) 


9? 


XRROrFLOAT ( IRRO) 


93 

9554 

CALL PLOKXRRO* Y,1 , I NO • 1 8 0 • 0 « - 160 . 0 ) 


94 


pho=?oo. 


95 


WN0RW=2,*CLXP(CMPLX(0, . -?.*PI*RHO/Wr ) )/SQRT(RHO ) 


96 


WRITEC6.311 ) 


97 

311 

FORMAT (1HI*7X,16H RELATIVE H FIEL0*33X.13hT0TAL H 

FIELD /2X « 9HMAGNI 

96 


?TU0E*6X t 5HANGI.E • 7 X t 2H0R , DX . 5HTHE T A , 9X , PHDB , 6 X . 

9HMAGNlTU(JE*7X. 

99 


3 5H ANGLE ) 


mo 


no 317 JNX=1»180 


101 


ThS=0 *017 4 5329 ♦FLOAT ( JNX ) 


102 


T=CMPLX(O.0*0.n) 


103 


no 310 I=1«M 


104 


yn=o.5oo*(x( i >+x< i+i) ) 


105 


HM=0.500*(H( I )+H( 1+1 ) ) 


106 


THN=TN ( I ) 


107 

310 

T=T+ ( (F( I )*CEXP( C v PLX(0.0<G*( (XN*COS(THS) ) *M 

HM*S I N ( THS ) ) ) ) ) ) 

106 


? »COS(THfj-THS) )»SEGL( I ) 


109 

r 

********** this CORRECTS THE OUTPUT TO TRUE MAG • 

FIELD 

110 


t=t*sts 


111 


rM=CAUS(T) 


IIP 


DB=2O.O*AL0G10<CM) 


113 


rANG = 57.P96*ATAN2( AIMAG(T) ,RFAl<T> ) 


114 


TMSD=TmS*57.296 


115 


abfs< JMX)=CM 


116 


TALL WlLLUM(20n. tTHSO*HlNCtWF) 


117 


HTC=T+HINC/HN0RM 


116 


WTM=CA°S(HTC) 


119 


HT A=RDEC* AT AN2 ( A I 'TAG ( HTC ) t REAL ( HTC ) ) 


120 


HTD=PO.tALOGlO(HTM) 


121 


wTOT< JNX)=HTM 


1 ?2 

317 

WRITE (6,312) CM f CANGtOB,THSDtHTU,HTMf HTA 


123 

312 

FORMAT(Fl 5 ,fN 2 ri 0 . 4 . 2 X«F 10 . 4 . 2 X,F 10 . 4 «L 15 . 8 tF 10 . 

4) 

124 


00 9500 JC=lil80 


125 


YO )=ARfS( JC) 


126 


Y ( 2 ) =HTOT ( JC ) 


127 


U=FLOAT( JC) 


128 


TNOrjC-J 


129 

9500 

CALL PLOT (U«Y*2tINDtl*OOfO«0) 


130 


PLOSS=O.U 
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131 

13? 

133 

134 

135 

136 

137 
136 

139 

140 


WRIT! («>« c » f ’9^) 

9994 ruWMAlC* Klr'fR LOSS f‘fR UNIT RtS ] ST I V 1 Ty * / • ft • « 5X * * PLOSS f ) 

nu 909*1 L USSrj f N 
TtMP.rnnsu (loss) ) 

PL0SS = )U.0SS4 Tfc MP*TFMr*SrGL(LOSS) 

WRI Tt ((it 9996 I i OSS * PlOSS 
9996 P0nMAT(I5 f F35.6> 

9999 ruNTUJUF 
TALL EXIT 
FNO 
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SUBROUTINE GDPLDV 


Subprogram GDPLDV is the groundplane surface division subroutine 
which is designed to symmetrically place 8 short segments around the source 
pair and divide the remaining segments along the rest of the groundplane. 

In use, this subprogram is called only once with the calling parameter NGP 
set to the number of matching segments one desires on the groundplane side 
of the origin. The 8 segments around the sources each have lenqth equal to 
0.25 x the distance between the sources. These segments are placed two to 
the left of the left source, four between the sources, and two to the right 
of the right source. The remaining NGP-8 segments are placed with at least 
one segment between the left endpoint and the symmetrically placed segments 
and the rest of the segments between the right of the symmetrically placed 
segments and the origin. The numbering of the endpoints is from the left 
endpoint and is indexed by variable MP. The cooridinate of the points are 
stored in [X(MP) ,H(MP)] and the surface normals and segment lengths in TN 
(MP) and SEGL(MP) respectively. These variables are shared with the main 
program and other subroutines through COMMON/MAINDV/. The subroutine listing 
follows . 
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: 

? 

3 

4 

5 

6 
7 
6 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 
3n 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 


subroutine GOPLnv(>ir,f>,rp) 

rOWMOi./MAlM/V/ X* ?3 0) ,Ht ?30) ,Tf„ ( ^30) ,SEGL (230) * DEL ♦ HC , 6 * CPL < EPR ' r,‘L 
common /r-AwiGO/ D2i,r>^2 
OtLl = APSlEPL/<MPP-4 ) ) 

Dtl2rA0S< ( DZ2-DZ 1 >/2.0> 

Xl=0z2-0( L2 
X2=022+0EL2 
X3=0?i-0tL2 
X4=021+0LL2 

WHITE <6, 20) OELl , 0EL2 * XI , X2 , X3, X4 

20 rORMAT(6H2.6> 

TU = ARS(EF J L~X1) 

MlsiriXCTlJ/OPLl) 

1MN1> 1 « 1 * 2 

2 nU.=TU/Nl 
no 3 1=1, M 
*m=rpL-Mi-n*OEi 

StGLt I >=UEL 

3 MP = I 

GO TO 4 
1 MP=1 

XlMP)=rPL 
SE.GL ( MP ) = TD 

4 UJ=MP 
DEL=0E:L2/2.0 
ro 5 1 = 1,8 

JSJJ+I 

y ( j)=xi+ < i-i ) *r>CL 
SEGL ( J ) =C»EL 

5 WP=J 

j JrMp 

N1=NGP-NP 
TF(Nl) 6,6,7 
7 ntL=ABS(X4)/M 
00 9 1=1, Ml 
J=JJ+I 

y ( J)=X4+( i-i ) *n£L 
SEGL<J)=DEL 
9 MP=J 
0 00 10 1=1, MP 
TNC 11=1.57079 
10 W < I ) = 0 . 0 

return 

6 WKITE<6,?1) 

21 PORMaT < • POOR POINT SELECTION* ) 

GO TO a 

rN o 
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SUBROUTINE VSURDV 


The subroutine VSURDV performs the surface division on arbitrary 
shape vee corrugations. Because of the periodicity of the corrugated sur- 
face, only one corrugation need be specified in detail and all the remain- 
ing corrugations are images. Figure A1 shows a period of the corrugated 
surface with the comer locations specified by the input parameters XI P to 
X5P . The input to this subroutine requires the corners be located in cm. 
rather than normalized to the corruqation period. Lines 28 to 32 of the 
main program remove this normalization. The segment endpoint numbering be- 
gins at "MP" at the origin and continues to the right endpoint. If the 
calling parameter KSURF is 1 (or 0) the surface will be closed (or open) 
and the numbering continued alonq the dotted path of Fiq. A2 (or the num- 
bering terminated). Some shapes which can be handled are shown in Fiq. A3 
and include the symmetric V-shape (case 1) and the square corruqation (case 
4 ) as well as asymetric and flat bottom vee corruqations . The subroutine 
1 i sting fol 1 ows . 

1 SUBROUTINE V/SUROV ( X 1 ♦ X2 * *3 , X4 , T C , NX , M2 , N3 • N4 . N5 • MP • KSURF ) 

PC IF K^URFr 1 , SURF ACT IS CLOSED 

3 C If KSURF=0, SUPFACF IS OPEN 

4 CQ^HON/MAINOV/ X(r>30) ,H(230) ,TN(230> fSECL (230) , OEL , HC , G , F PL * EPR , WE 

5 OATA Pl,Pl2/3. 1415926, 1.570796,5/ 

6 w = 

7 JU=MP 

B T F ( N 1 ) 2*2,1 

9 1 PEL=X1/N1 

in ' no in t = l,M 

11 j=JJ+I 

12 X(J)=DEL*FLOAT ( 1*1 ) 

13 

14 TN(J)=PI2 

15 SEGL(J)=OEL 

16 MP=J 

17 10 CONTINUE 

18 2 xl=X(«P)+OFL 

19 nX=CX2-Xl)/N2 

20 nH = HC/r>2 

21 nLLrSORT { OX *OX + OH * OH ) 

22 T=ATAN?(UX,DH) 

?3 jj=mp 

24 no 20 1=1, N2 

25 ,I = JJ+I 

26 V(U)=XT4DX^FLOAT( 1-1 ) 

27 H< J)=-DH*FLOAT ( 1-1 ) 

2ft TN(J)=T 

29 SEGL(J)=UEL 

30 MP= J 

31 20 CONTINUE 

32 VT=X(MP)+DX 

33 JJ = * V 1P 

34 IF ( N3 ) 4,4,3 

35 3 PEL=(X1-X2)/N3 

36 HO 30 1=1, N3 

37 J = JwHI 

36 X ( J)=XT+QEL*FLOAT( 1-1 ) 

39 M(J)=-HC 

40 TN<J)=PI/ 

41 RlGL(J)=DEL 

4 P rp rj 

43 30 CONTINUE 
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44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
6? 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 
74 

74 

75 

76 

77 
76 

79 

80 
81 
82 

83 

84 

85 

86 
87 
86 

89 

90 

91 

92 

93 

94 

95 

96 

97 
9fl 
99 

ino 

101 

102 

103 

104 

105 

106 
107 
106 

109 

110 


XT =X ( MP ) + DEL 

4 nx=(X4-X3)/N4 
nH = HC/\'4 

OLL^SQRT <DX*OX+DH*OH) 
T=1*57079*AT AN2tDH,0X) 
JJ=MP 

ro 40 I — 1 » N 4 
J=JJ+I 

X( J)=XT+DX*FLOAT(I-.l) 

W( J)=-HC+0H#FLPAT(I-1 ) 

TNC J)=T 
SEGL < J ) -DEL 
MP=J 

40 CONTINUE 

XT=X(FP)+OX 

JJ=MP 

TF(Nb) 6*6*5 

5 nLL = <TOX 4 )/N5 
no 50 I~1*N5 
J=JJ+I 

X(J»=XH DLL ♦FLOAT ( 1-1) 

W( J>=0*0 
TNI J)=PI2 
StGLC J)=OEU 
MP=J 

50 CONTINUE 

6 NTCELL=IFIX( <fpr+o«i >/tc> 
NCEUL=N1^N2+N3+N4+M5 

no 60 N=2,NTCELL 
JJ=M+NCFLL*<N-1) 

00 61 I rl * NCELL 

J=JJ+I 

K=M + I 

X(J)=X(K)+TC*(N-1) 

W| J)=H(K) 

TN < J ) =TN< K ) 

SEGLC J>=SEGL<K) 

61 ^P=J 

60 CONTINUE 

TF(KSURF) 7*7*8 

7 MPrP'JP + l 
X<NP)=NTCELL*TC 
u(MP)=0.0 
UKlTF_(6*t,2) 

62 ruRM AT ( / ' OPfN SURFACE*/) 
return 

8 XMAXrNTU LL*TC 

.IJ = MP 

JJ = «P 
PP=HC+0.1 
M=IFIX(HP*8./WE) 
oel=hp/n 
no 70 I=1*N 
J-JJ+ I 
XI J)=X|-1AX 
M(J)=-(I-l)*OEL 

TN< J)=n.O 
c£GL< J)=UEL 
70 MP=J 
JJ=*IP 

K=IFIX( (XMAy-EPL)*8./WF) 

nEL=(XHAX-EPL)/K 

DO 80 I=1*K 

x(J)=X M AX-( I-l)*OEL 
W( J)=-MP 
TN( J)=-PI2 
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I 1 I 
1 1 ,* 

I I ^ 
1 )*» 
If* 
111. 
117 
116 

119 

120 
121 
122 

123 

124 

125 

126 
127 
126 

129 

130 

131 


st r.if J)sh» L 
an M'-j 
ij mi* 

r - h j x ( mi * t . /ur > 

n f i -nl*/fi 
no 9 n J - ) # M 

J = vJJ* i 

x< j>=lpl 

p< jIs-hpmi-i )* mu 

TN< J)=P1 
StGU J)=Ut L 
90 MP=J 

X <MP)=EPl 
H("P>=0.0 
TN(MP)=O,0 
StGLtMP>-U,0 
WHlTE(6«9l ) 

91 FORMAT ( / f closed SURFACE*/) 
return 

FNn 
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SUBROUTINE CO 


Subroutine CO rjenerates the coupling coefficients for the "v H 
matrix of section II. The subroutine finds the magnetic field at seqment 
"MR" duo to a current on seqment "MC" by performing a 5-point Gaussian in 
tegration along the length of segement "MC". The program listing follows 


i 

? r 

3 

4 

5 

6 
7 
6 
9 

10 

11 

12 

13 100 

14 

15 

16 

17 

18 
19 
?0 
21 
22 

23 

24 

25 

26 
27 
2 A 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 


complex function CO < MR « Mf ) 

THIS GIVES THE OLD MATRIX COEFFICIENTS 

rOMMON/MAINDV/ X ( 230) ,H(?30 ) ,Trg(?30 ) ,$EGL (230) • DEL • HC * 6 • EPL « EPR i UE 
COMPLEX UJC 
COMMON /COG/ DJC 
COMPLEX AHAN21 

DATA Gul • Gu2 •Gu3tG U 4 i GuStGWl *Grf2»GW3tGw4 ,Gw5/-ft. 9061798 t -0.5364693 

2,0.00,0.5384 693, 0,9061798, 0.2369266, 0.4 766287, 0.56 9868, 0.4786287, 

3 0.2369268/ 

TF(MR.PE.MC) GO TO 100 
f 0=CMPLX (0.500,0.0) 

CO TO 200 
CONTINUE 
DEL=SEGL(MC) 
nEL2=DEL/2 . 

XMM=0.50U*(X(MR)+X(MR+1 ) ) 

MXMM=0.50t (H(M^ )+H(MK+l ) ) 

FPL=X(MC) 

FPU=X(MC+1) 
nVOFFP=(tPU-FPt )/2.0 
nVSMEP=(t PU + FPL ) /2.0 
yU5=GU5*UVDFLP+0VSMEP 
yUl = GUl*DVL)FEP + OVS v EP 
yu2 = GUP*L , VL'FEP + f’VSMrP 
XU3=GU3*UVDFEP+DVSMEP 
yU4=GU4 ♦UlLFEP+OVSMEP 
HU=H(. w IC + 1> 

HL=H(MC) 

nELH= ( HU-HL )/2. 

nHM=(HU+HL)/2.0 

HXU1=GU1*0LLH«m>HM 

pXU2=GU2*ULLH+nHM 

HXU3=GU3^DFLH+nHM 

HXU4=GU4*0ELH+nHN 

MXU5=GU5*bELH+nHH 

THN”TN { MC ) 

XR=XMM-XUl 

YR=HXMM-HXU1 

P10N=( XR*COS(ThN)+YR*SIN(ThN) ) /SORT ( XR* XR+YR*YR ) 

yR=XMM-XU2 

YR=HXMM-HXU2 

R2DN=(XR*C0S(TMN)^YR*SIN(THN) >/SQRT(XR*XR*YR*YR) 

XR=XMM-XU3 

YR=HXMI w -hXU3 

P30N=(XR*C0S(THN)+YR*SIN(THN) ) /SORT ( XR*XR+YR*YR ) 

yR=XMM-XU4 

YK=HXM m -HXU4 

R4MN=( XR*COS(TwN)+YR*$IN(THN> ) /SCRT < XK*XR+ YR *YR ) 

XR=X MM-XU5 
YR=HXMM-HXU5 

R50N=(XR*C0S( TMN)+YR*SIN(lHN) )/SGRT(XR*yR+YR*YR) 

PACT oR = SOK T (OV n FEP*OVDFEP + OEl H*OELH ) 
co=factor*( 

?+(GWl*AHAN21 (G*SQRT( ( ( XU1 - X*" ) **P ) ♦ < ( HXtJl -HX^w ) *♦? ) ) >*RlON ) 

?+{ GW? * AM AN? 1 (G*SQRT ( ( < XU2- Xwm )**?)♦ ( ( HXU2-HXMK ) * *2 ) ) ) *R?DN ) 
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56 

57 
5P 

59 

60 
61 
6? 


2-MGW3*AH£f'?l (6*Sf'RT < ( (XU3-XM M )**2) + < (HXij3-HXm*i)*% 2) ) ) ♦R^Df) 
?+(GWM*AHA'.?l(t*5fy^T C ( ( X(J4 -X*m ) **£ ) ♦ < (H> ‘ ) **2 l ) ) *F4U^) 

?+(GW5*AHAfJ?l(6*$( 1 iRT ( ( ( XU5-XP™ >**?)♦( (HXU5-HXF ''>**?> ) ) *R5DN ) ) 
rO = C r >* f lJC 
200 rOMTlNHt 
RETURN 
TND 


45 


SUBROUTINE AHAN21 


Subroutine AHAN21 is a double precision Hankel function of type 
2 and order 1. The function is generated using the polynomial approxi- 
mations for J] (X) and Yi(X) presented in Abramowi tz[13]. The program 
listing follows . 


1 

2 r 

3 

4 

5 

6 
7 
b 
9 

10 

11 

12 

13 

14 

15 

16 
17 
10 

19 

20 
21 
22 

23 

24 

25 200 
26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 300 

42 
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FUNCTION AHA N 21 (X) 

THIS IS THE HANKEL FUNCTION OF TYPE 2 A NO OP ORDER 1 
DOUBLE PRE.CISIO 1 '' X0tDX«Al,A?, A3,A4,A5,At,HJl f ^l»R2 f B3 f R4tB5, AHJ1, 
2T0X, Al # A? ( A3,AU , A5 % A6 f Tl f T2 f T3 , T4 , T5 , T6 , T7 , UGQX , "6 
COMPLEX AHAM21 
nX=OOLP (X) 

IF IX.GT.3.0) GO TO 200 

yu=ox*nx/°. oo+oo 

Al = - 0 . 3176 in- 0 X + 0 .nU 9 P -0 4 *xn 

A2 = C. 0044331 90+00+ A1*X0 

A3=-0*n 39542 89P + 0 0-t A2*XD 

A4=0 .21 093573P+00+A5*vD 

Ab=-0.56249965n+00+A4*xn 

Ab=0.5n+00+A5*Xn 

MJ1=A6*0> 

Pl=- 0 . 040 097 60 + 0 0 + 0 . 0 027873 D+ 0 0 *X 0 

P 2 = 0 . 31239 MO + 0 0 + 31 *XO 

P3=-l . 3l64A270+00+n2*xn 

P4-2 *16827090+ 00 + 93* XH 

*5=0. 22 12091 0+00+94* XU 

P6=- 0.63661 9 60+n0+«5*xn 

AHJ1 = (P6/Dx) +H.M 4 OLOG( OX/2.0) *0*63661977 
AHAN21 =CMPLX (S m GL(HJ 1 ) ,-SNGLf ANJ1> ) 

GO TO 300 
TUX= 3 . 0 /DX 

A 1 = 0 . 0011 3653 D+ 00 - 0 . 00020033 *TUX 

fl2=-0. 0024951 1H4 CO +A1+TPX 

A3 = .000171 O5n+O0+A?*TOX 

A 4= 0.01 65 96670+00+ A 3* TPX 

A5=0. 1 560-05+ AU+TOX 

Afa=0.7P7P 545 6^+00 +A5*TDX 

Tl=O.OrtO 79624 0 +O 0 - 0 . 0 n 0291 £ 6 P+ 00 *TOX 

T2=0.00074348P+on+Tl+Tnx 

To=-0.00637C79 n +00+T2*TOX 

T4=0. nnOU66bOO+OC+T3*TPX 

Tb=0. 1249963 20+0 0+T4*TPX 

T 6 =- 2 . 3561944 90 + 0 0 +T 5 *TQX 

T 7 = 0 X+T 6 

PSQX= Afa/DSCRT (nX) 

AHAN 21 ^CNPLX(S^'GL(DSQX* 0 C 0 S(T 7 ) ) t-SNGL(OSQX* 0 SIN(T 7 ) ) ) 

CONTINUE 

return 

end 
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SUBROUTINE CROUTN 


The nonsymmetric matrix inversion subroutine solves a system of 
simultaneous linear equations with complex coefficients using the method 
developed by P.D. Crout[10]. The set of equations are related to the call 
ing parameters shown in the following listing by 

[C] [X] = [F] 

where [C] is the array of complex coupling coefficients 

[X] is the column vector of length N of unknowns, 

[F] is the column vector of length N of known variables before in 

inversion; after inversion, the array contains the solution 
vector. 

NDIM is the dimension of [C] and [F] in the subroutines; NDIM >.N. 
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1 SUBROUTINE. CROUTN ( C » F , M « MuIM) 

2 COMPLEX C ( NOI M , NO I M ) » F ( NO I M ) iS 

3 C NONSYMMF.TRIC CROUT 

4 C FIRST COLUMN OK 

5 c to get the first row 

6 no 10 J = 2.N 

7 10 r ( i,j)sC(3 ,J)/r<i,u 

ft r NOW WORK ON ROW AflO COLUMN SET K 


9 no 11 K = 2,N 

10 KM0=K-1 

11 KP0=K+1 

12 C TO GET DIAGONAL ELEMENT 

13 SsCMPLXIO, 0,0.0) 

14 NO 1? IK=J , KMO 

15 12 S=S+C(K,IK)*C<TK,K) 

16 C (K,K)=C(K,K)-S 

17 C TO GCT ELEMENTS IN COLUMN K BELOW ROW K 

16 TF (KPO.GT.N) GO TO 17 

19 DO 13 IROWsKPO.N 

20 S=CMPLX(0. 0,0.0) 

21 no 14 JJ=1,KM0 

22 14 S=S+CUR0W«JJ)*C( Jd.K) 

23 13 r»IROW.K)=C(IROW,K)-S 

24 C TO GET ELEMENTS IN ROW K TO THE RIGHT Of COLUMN K 

25 no 15 ICOL=KPO,N 

26 S=CMPLX(0. 0,0.0) 

27 no 16 JR= 1 ,KM 0 

26 16 S=S+C(K,JR)*C<JR«ICOL> 

29 15 C(K,ICOL)=(C(K,ICOL)-S)/C(K,K) 

30 17 CONTINUE 

31 11 CONTINUE 

32 C THIS ENOS THE MATRIX FACTORIZATION 

33 C THIS BEGINS THE BACK SUBSTUTION 

34 C CONVERSION OF SOURCE SIDE 

35 F(1)=F<1 1/CI1.1 ) 

36 no 90 IJ=2,f* 

37 SsCMPLX (0.0, 0.0) 

30 TJM0=IJ-1 

39 no 91 IK = 1 ,IJM 0 > 

40 91 S=S*C( IJ, IK)*F ( IK) 

41 90 F< IJ)=(F( IJ)-S)/C( IJ.IJ) 

4? C NOW FOR FINAL BACK SUBSTITlTION 

43 MM0=N-1 

44 00 160 L=1«NM0 

45 K=N-L 

46 KP0=K+1 

47 S=CMPLX(0. 0,0.0) 

46 HO 175 jOsKPO , N 

49 175 S=S*C(K,J0)*F(J0) 

50 160 F(K)=F(K)-S 

51 C THIS ENDS THE BACK SUBSTITUTIONS 

52 RETURN 

53 END 
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SUBROUTINE WILLUM 


Subroutine WILLUM is the subprogram used to generate the incident 
fields at RflO, Pill from a magnetic line source pair located at U/l and L)/?. 
The line source amplitudes are 1/Z 0 and phased to produce a cardiod pattern 
directed along 0=0 (toward the corrugations). The Fortran listing follows. 

1 SUBROUTINE WILLUM! RHOiPHI «himc*wet ) 

2 COMMON /H N AV‘IGl>/ UZi t DZ2 

a r this subroutine, calculates the near fields at rho.phi Radiateo by a pair 

4 r WIRES SEPARATED BY o AMO H ABOVE ORIGIN WITH MAGNETIC CURRENTS FK 1 ANo FK 

5 COMPLEX HlNCtFKl«FK2t JN ( 2 * 6 ) « J 

6 COMPLEX aRG 

7 PI=3. 14159 

8 nEGR = PI/16(j. 

9 FK=2.*PI/SE 

10 FKl=(l.tO.) 

11 JsCO.tl.) 

12 FK2=J 

13 H=0. 

14 DZl=-4. 

15 nZ2=-6. 

18 RH01=SQRT (DZ1*DZ1+H*H) 

17 RH02=S0RT ( 07? ♦ OZ2+H*H ) 

18 PH1=ATAn?(H.OZ1 > 

19 pH?=ATAN?(Ht072> 

20 RHOlPrSOPT ( APS ( RHO»KHO+RH01*RH01-2. ♦RHO*RHOl*COS ( PHl.PHl *DEGR ) ) ) 

21 PH02P = SQRT I APS I PHO *Kh0+PHO2*RH02-2 . *RhO*PHO?#cOS ( PH2-PHI*DEGR ) ) ) 

22 AKG = CMPlYIP K»RH01P#0* ) 

23 IF (CABS(ARG) •GF.O.Ol) GO TO 1 

24 H1NC=FK1 

25 no TO 2 

28 1 CALL BESSEL( ARGt JN) 

27 hIMC=F^1»JN(1,S) 

28 ? CONTINUE 

29 ARG=CMPLY ( FK*RH02P , 0. ) 

30 IF ( CABS ( ARG ) *Gf • 0 • 01 ) GO TO 3 

31 HlNCrFK2 

32 GO TO 4 

33 3 CALL BESSEL! ARG* JN) 

34 MlNC=HUJC + FK2i'JN(l ,5) 

35 4 CONTINUE 

38 C 20=377. SCALE FACTOR REMOVED 

37 HlNC=HlNC*FK/4, 

38 RETURN 

39 FND 
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SUBROUTINE BESSEL 

RESSEL generates J 0 (x)> N 0 (x), Ji(x), Ni(x), Ho^x), h| 2 |x) and 
their derivatives using a power series for small complex arguments, |7.| c 12, 
and the asymptotic forms for |7.| > 12. These variables are stored in array 
JN and returned as a calling parameter. The Fortran listinn follows. 

1 SUBROUTINE PESSELI ?, JN ) 

2 (-(IMPLEX JN<2.6) .Z. TERMJ,TERMM.MZ24.T1 .T2.T3 

3 r JN<1,1»=J0(7) 

4 r / 

5 r jn<i ,<j)=m)(Z) 

6 C JN(2.2)=Nl<Z. > 

7 r jfjci,a»=i)/r-x<j0(xn 
6 c JN« 2 ,a)=ivnx(Ji(xn 

9 C JlM(l,t»)=0/»'X(^i)(X) ) 

10 c JN(?. 4 )=n/nx(Ni (x> ) 

11 C JN<l,b)=H0<2) (X) 

1? r JN (2,5) —Hi (2 > ( X > 

1? C JN(l.b)=O/PX(H0(2» ( X » ) 

14 r JN< 2 , 0 >= 0 /OX(H 1 (2) (X) ) 

15 PI = 3. 1411)9 

IS TF(CAHS( z ) . GF • 1 2 . 0 ) GO T010 

17 rACTOR=n.O 

in TtRM(g=(n, ,o. ) 

1 9 MZ24--0.2 L*?*Z 

20 tERMj- ( 1 , 0 , 0 .0 ) 

21 nu 1 Nl*=1.2 

22 N=MP-1 

23 jNJNP.l > = TERMJ 

24 M = 0 

2b 2 m=M+1 

26 TERMJ=TEKMJ*MZ24/FL0AT(M*(N+M) ) 

27 ,|N(NP» 1 )=JN(NP, 1 ) + TF.RMJ 

26 IF(NP.NF.l) GO TO 3 

29 FACTOR=FACTOR+1.0/FLOAT ( M ) 

30 TtRMN=TERMM+TERMJ*FACTOR 

31 3 FRKOR=CABS( TERMJ) 

32 IF (ERROR, GT.l.OE-lO) GO TO 2 

33 1 TLRMJ=0.S*Z 

34 JN 1 1,2 ) = <2. 0/3.14] 5927 )*( C 0 . 57721 5?+CLO(, < 0 . 5*Z ) )*JN<1 ,1 ) -TERMN ) 

35 11 rONTINUE 

36 JN(2,2>=( JN(2»1)*JN(1,2) -2.0/(3.1415927*2) 1/ JN(1.11 

37 JN < 1 , 3 ) =■ JN ( 2 , 1 ) 

36 jN<2.3)=JN<1.1 >-JM2.1)/Z 

39 JN (1.4) = - JN (2.2) 

40 .IN (2.4 I = JN (1.2 t * JN < 2 > 2 ) /Z 

41 JN( 1 . 5 ) =JN (1.1 )-(0.0,1.0)*jN(l,2) 

4? JN ( 2 , 5 ) = JN ( 2 » 1 )-( 0. 0,1.0 )*JM( 2.2) 

43 JN (1,6) = JN ( 1 . 3 ) * ( 0.0,1 .0)*JNd .4) 

44 jn( 2,6)=JN(?,3»-( n.0,1.0)*JN)2,4) 

45 RETURN 

46 c asymptotic forms usfd for arguments greater than 12 . 

47 10 Tl= Z -CMPLX(Pl/4.,0.) 

46 T tt- Z -CMPLX(3.*PI/'4.,0.) 

49 T3=CS0RT(2./(PT* 7 ) ) 

50 .IN(l.l)=Ti*CCOS(Tl ) 

51 ,)N(2.1)=T3*CrOS<T2) 

5? JN(1.2)=T3*CSIM(T1 ) 

53 JN(2,2)=T3*CSIM(T2 ) 

54 GO TO 11 

55 ENO 
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SUBROUTINE PLOT 


PLOT is a line printer plotting subroutine which plots up to 10 
different curves (stored in the Y array) on the same set of coordinates. 
The value "X" is plotted along the page while the "Y" array is plotted 
between YM1N and YMAX across a line on the page with Y(l) plotted as a 
Y(2) plotted as etc. A call with "IND" = 0 places the Y scale 
across the page and for IND equal to a multiple of 10, places tic marks 
across the page. The PLOT subroutine must be called in a "DO" loop to 
plot an array of data. The Fortran listing follows. 

1 4URROUT I [jfc.Pl OT ( X«Y.N«mt»YMAX,YMIN) 

? IYI PENS loop ( 1191 » Yt/'PEL(6) . Y(10) .MARK <10 ) 

3 PAT A MARK ( 1 ) ,MARK<?> .MARK (3) , MARK ( b ) • MACK ( 6 ) , MARK ( 7 > , MARK ( 8 ) « 

4 2 MARK(9» , MARK (10) ,MARK(4)/1H*,1H. ,1HJ,1H0,1HN,1HH,1H1,1HZ,1H-,1HX/ 

5 nATA IBLAMK.NOPT,IPlUS/1H « 1H* « 1H+/ 

6 TF < INO) 1 . 1 , 11 

7 1 W R I T E { 6 , 3 ) 

6 3 f0RMAT(lHl//25X*4BM0R0FR IN WHICH PLOT SYMBOLS APE USED *.IXONhlZ 

9 ♦ <-//30X.39hTHE SYMBOL (*) INDICATES OFF-SCALE DATA//) 


10 


D07J=9,119 

11 

7 

i*M J>=MARK ( 10 > 

12 


MCOUNT = lfJ 

13 


«CALE=lO0. 0/( YMAX- YM IN) 

14 


1 ( LL=(-YM1N*SCALE>*1 1.5 

15 


noej=i *b 

16 


p=J-i 

17 

A 

YUABF.U J)=R*20.0/SCALE + YMIN 

1 A 


WRITE (6, 9) < YLARFU I ) * 1 = 1. 6) 

19 

9 

f-ORMAT (6X , 1 P F 9 * 2 t 5 ( 1PL20.2) / ) 

20 


6 0 T 0 1 3 ? 

21 

11 

NCOUNT = NC OUNT + 1 

22 


D099j=3 ,119 

23 

99 

M( J) sIRLANK 

24 


tMLLL.GL.il .ANO.LLL.LE. 110) NULL )=MARK (10) 

25 


tF(NCOUNT-lO) 13?.l32i 133 

26 

13? 

D069J = 11 • 1 1 1 « 20 

27 

69 

M( J)rIRLUS 

26 

133 

0020j=liN 

29 


tsC Y( J)-YNIN)*5CALE+0.b 

30 


fF(L)14, 17.17 

31 

14 

T F <L+10)15, 16,16 

32 

15 

Ml) sNOPT 

33 


6OTO20 

3*4 

1 6 

LL=L+11 

35 


M(LL)=HARK( J) 

36 


G0T020 

37 

17 

tF (L-1O0)16,19,19 

36 

16 

LL=L+11 

39 


M(LL)=MARK( J) 

40 


gotopo 

41 

19 

M119)=N0PT 

42 

20 

CONTINUE 

43 


t F ( NcOUNT - 1 0 )21.25,21 

44 

21 

WRITF<6,?4) (K,(J ), Jr 1 , 1 19 ) 

45 

24 

rORMAT(lX,119Al) 

46 


G 0 T 0 ? 7 

47- 

25 

ViHlTf (6,26) (X, tut J) • J=9.1l9> ) 

46 

26 

t ORM A T ( 1 X , F6 • 3 « 1 1 1 A 1 ) 

49 


NCOUNT = n 

50 

27 

CONTINUE 

51 


RlTURN 

52 


F*ND 


51 


APPENDIX B 


HIGHER ORDER MODES IN A SQUARE CORRUGATION 

Because of the limited number of matching points available, only 
four segments were used on each of the vertical walls of each of the cor- 
rugations. Over the range of corrugation depths considered, this meant 
that 8 to 16 matching points per wavelennth were used. While this den- 
sity is adequate to represent the TFM mode in the square corruqation, any 
appreciable amplitude higher order mode would either not be observed or 
could introduce errors in the computations. Thus it was necessary to verify 
that any higher order modes were of small amplitude and only over a limited 
portion of the corrugation. To do this, the closed surface of Fia. B1 was 
used with the computer program of Appendix A. The model included a sinqle 
corruqation illuminated by a magnetic line source located n.5x (4 cm) from 
the corrugation. 

The 30 segments per vertical wall in the corrugation yields 60 to 
240 matching points per wavelenqth over the ranqe of depths considered 
(0.125 d/x < 0.5). The relative surface current which exis.ts at each 
point on the surface of the corruqation is shown on the plots in Fiq. B2a 
to B2d for corruaation depths of 0.125X, 0.25x, 0.375X and 0.4Qn\ respectively. 
On these plots, the surface current amplitude, normalized to the l/7 0 amp- 
litude magnetic line source, is plotted along the corruqation wall where 
the current is flowing. The amplitude function cos ( d-H) / x corresponding 

to that of the TEM mode whose maximum is set equal to the maximum current com- 
puted in the corruaation is also shown. Notice that in each case the calcu- 
lated current (the series of steps) very nearly fits the cosinusoidal distri- 
bution anticipated for the TEM mode in the corrugation. Also, the computed 
currents on the two walls at the same depth differ only slightly. Thus the 
higher order modes are of low amplitude amd only over a small reaion of the 
corrugation. The loss associated with higher order modes is negligible. 
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H = —4cm 34 SEGMENTS 


Fig. B1 . Surface model used to study hiqher order modes in a 
single square corrugation. 



Fig. B2a. Relative surface current existing on the walls of a 

square corrugation at a corrugation depth of d/A = 0.125. 
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20 SEGMENTS 



Fig. B2b. Relative surface current existing on the walls of 
a square corruqation at a corrugation depth of 
d/A = 0 . 250 . 
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Fig. B2c. Relative surface current existing on the walls of a 

square corrugation at a corrugation depth of d/x = 0.375. 
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Fig. B2d. Relative surface current existing on the walls of a square 
corrugation at a corrugation depth of d/\ = 0.499. 
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